The s — | Antiferromagnetic Heisenberg Model on Fullerene-Type Symmetry Clusters 
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The Si = | nearest neighbor antiferromagnetic Heisenberg model is considered for spins sitting on 
the vertices of clusters with the connectivity of fullerene molecules and a number of sites n ranging 
from 24 to 32. Using the permutational and spin inversion symmetries of the Hamiltonian the low 
energy spectrum is calculated for all the irreducible representations of the symmetry group of each 
cluster. Frustration and connectivity result in non-trivial low energy properties, with the lowest 
excited states being singlets except for n = 28. Same hexagon and same pentagon correlations 
are the most effective in the minimization of the energy, with the n = 32 — D^u symmetry cluster 
having an unusually strong singlet intra-pentagon correlation. The magnetization in a field shows 
no discontinuities unlike the icosahedral //, fullerene clusters, but only plateaux with the most 
pronounced for n = 28. The spatial symmetry as well as the connectivity of the clusters appear to 
be important for the determination of their magnetic properties. 

PACS numbers: 75.10.Jm Quantized Spin Models, 75.50.Ee Antiferromagnetics, 75.50.Xx Molecular Magnets 



I. INTRODUCTION 

The antiferromagnetic Heisenberg model (AHM) has 
been the object of intense investigation for some time 
now as a prototype of strongly correlated electronic be- 
havior. The effects of frustration, quantum fluctuations 
and low dimensionality can lead to new phases differing 
from conventional order and possessing a non-trivial low 
energy spectrum 0, 0, 0| . Small magnetic clusters pro- 
vide an excellent testing ground for the validity of the 
AHM as well as other theoretical models, as oftentimes 
its low energy properties can be computed on these struc- 
tures and itspredictions can be directly tested against 
experiments [J]. 

The fullerenes are a class of three-fold coordinated 
molecules consisting of § — 10 hexagons and 12 pentagons, 
with n their number of vertices (or sites) Q. With in- 
creasing n their shape resembles more and more the hon- 
eycomb lattice, albeit in closed form, with the pentagons 
playing the role of structural impurities. Frustration de- 
creases with n, as the unfrustrated hexagons dominate 
in number the frustrated pentagons, while the distribu- 
tion of the latter determines the symmetry group of the 
molecule. The properties of the AHM have been com- 
puted for the smallest element of the family, the dodeca- 
hedron, which consists only of pentagons and belongs to 
the icosahedral point symmetry group Ih- For classical 
spins the signature of frustration is very strong, gener- 
ating three magnetization discontinuitites in an external 
magnetic field, unexpectedly for a model lacking mag- 
netic anisotropy [1, 0] ■ In the full quantum limit where 
the individual spin magnitude s, = i, the low energy 
spectrum consists of singlets, absent in unfrustrated sys- 
tems |8j. More unconventional behavior is displayed by 
the magnetization which is discontinuous in an exter- 
nal field as in the classical case, and the specific heat 
which has a two-peak structure as a function of tem- 
perature. For Si = 1 non-magnetic excitations are still 



present inside the singlet-triplet gap, and now there are 
two magnetization discontinuities in a field. Similar be- 
havior with magnetization discontinuities was also found 
for the AHM on larger fullerene molecules of Ih symme- 
try for classical and Sj = h spins [7|. There is strong 
evidence that this behavior is persistent for J/j symmetry 
in the n — > oo limit and survives asymptotically close to 
the zero and saturation fields, even though the number 
of hexagons strongly dominates the 12 frustrated pen- 
tagons. 

Motivated by the non-trivial spectral and magnetic 
properties of the AHM on the //, clusters, the Sj = \ 
AHM on fullerene clusters with symmetry other than Ih 
is investigated. Modine and Kaxiras calculated ground 
state energies and correlation functions of clusters with 
n up to 28, and by truncating the Hilbert space of a 
cluster with n = 32 [9!]. They used no symmetries to 
reduce the computational requirements and found that 
same hexagon correlations are the most important for 
the minimization of the energy. For the n = 28 clus- 
ter more specifically the ground state was found to be 
doubly-degenerate. Their results are extended here with 
the calculation of the low-energy excitations and the re- 
sponse in a magnetic field for the clusters of Ref. 
with n = 24 and 28, the most symmetric for these n 
values. The magnetic behavior of different clusters with 
n = 26 and 32 is investigated, as well as for a cluster with 
n = 30. All of them are shown in Fig. Q] [a.flb|. They are 
the most symmetric isomers for each n (for n — 32 the 
two clusters belong to different point groups who have 
the same number of symmetry operations). The spatial 
symmetry group along with the number of symmetry op- 
erations is listed in Table U for each n. The low energy 
spectra are calculated with Lanczos diagonalization, and 
point-group and spin-inversion symmetries are used to 
reduce the computational requirements and to classify 
the states according to the total symmetry group's irre- 
ducible representations 0] . Similarly to the /^-symmetry 
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case it is of interest to look for non-magnetic states in- 
side the singlet-triplet gap, and search for unconventional 
behavior of the magnetization and the possible presence 
of discontinuities in a field. Comparing with the case of 
the Ih clusters we gauge the effect of symmetry on the 
behavior of the AHM on the fullerenes. It is noted that 
fullerene clusters are edge- and not corner-sharing, and 
it is not obvious if it is possible to re-write the Hamilto- 
nian as a sum of total spins on individual units or how 
to perform any other mathematical operations in order 
to derive analytic results, even in the classical limit. 

From now on sites which belong only to pentagons will 
be called pentagon sites, while the rest will be called 
hexagon sites. The clusters will be abbreviated as C„, 
with the D^d n — 32 cluster C32,/, and the D$h n = 32 
cluster Cz2,n- Looking at Fig. [TJ the pentagons form a 
band in the middle of C24, while the hexagons a band 
in the middle of C30 and C32,/. For C26 and C28 no 
hexagons are sharing edges, while for C32.11 the hexagons 
are adjacent to each other in two groups of three. Close 
proximity of polygons of the same kind can in principle 
minimize frustration, which also decreases on the aver- 
age with n. As seen in Table U and Fig. [2] the ground 
state energy per spin closely follows the aforementioned 
arguments. The clusters with n = 30 and 32 have the 
lowest energy per spin, as the hexagons approach each 
other more and more with increasing n. C24 has the next 
lowest energy with the band of pentagons separating the 
two hexagons. For C26 and C28 with the disconnected 
hexagons the energy is maximum, even higher than the 
dodecahedron (C20) energy per spin, which is also plotted 
for comparison [8j]. 

Frustration also results in non-magnetic excitations in- 
side the singlet-triplet gap. The first excited state in all 
cases except C28 is a singly or doubly degenerate sin- 
glet (Tables HEIIl Fig. E]). For C 30 and C 32 j the sec- 
ond excited state is also a singlet, non-degenerate. In 
contrast, C2s's first excited state is a triply degenerate 
triplet. For C24, C30 and C32,/ the low-lying singlets 
are well separated from the magnetic excitations. For all 
the clusters the dominant nearest-neighbor correlations 
in the ground states are of the intra-hexagon type, with 
intra-pentagon following closely (Table IVIIip . The only 
exception is C32,// where the intra-pentagon correlation 
is very strong and equal to -0.60548. This is appointed to 
the geometry of the cluster, with the "strong" bonds lit- 
tle frustrated in the middle of the structure and between 
the two groups of three hexagons each. 

The magnetization in an external magnetic field does 
not present a discontinuity, as is the case for the do- 
decahedron Q. There are only magnetization plateaux, 
where the ground state belongs to specific total spin S 
sectors for a wider range of fields than their neighboring 
sectors (Fig. @|. The most pronounced case is for n = 28 
and S = 2. It appears that reducing the symmetry from 
Ih works against magnetization discontinuities. 

The clusters considered in this paper have relatively 
small n and mostly belong to different spatial symmetry 



groups. Only C26 and C32JI have the same spatial sym- 
metry, neverthless they have no common pattern of low 
energy and magnetic behavior as was the case for the Ih 
symmetry clusters 0, H|. Similarities are only found in 
the low energy behavior of C30 and 6*32,7 even though 
they belong to different symmetry groups, and they only 
have in common all the hexagons forming a band in the 
middle. It is also pointed out that there are some minor 
differences between the type I and type II molecules as 
they were called in Ref. Q , even though they both share 
the icosahedral Ih symmetry. Larger clusters can shed 
light on the magnetic properties as a function of the dis- 
tribution of the pentagons, however this is not possible 
with present day computational means, at least for the 
low energy properties. Neverthless, connectivity appears 
to be as important as spatial symmetry for the magnetic 
properties. 

The plan of this paper is as follows: in section [IT] the 
model and the method are introduced, and in section IIIII 
the low energy spectra and nearest-neighbor correlation 
functions are presented. Section [TVl presents the results 
on the ground state magnetization, and section [V] the 
conclusions. 



II. MODEL AND METHOD 

The antiferromagnetic Heisenberg Hamiltonian with 
spins Si located on cluster vertices i is 



H 



J E 



hS z 



(1) 



<»J> 



where <> denotes nearest neighbors, and J is positive 
and will be set equal to 1 from now on, defining the unit 
of energy, h is the strength of an external magnetic field 
and S z the projection of the total spin along the field 
direction z. For an explanation of the use of symme- 
try operations to block-diagonalize the Hamiltonian, see 
Ref. [8| and references therein. The data for the symme- 
try groups was taken from Ref. 11]. Degeneracies are 



reported with respect to states with specific S, each of 
which corresponds to 25 + 1 states with different value 
of S z . Lanczos diagonalization was performed in double 
precision, but results are shown with a smaller number 
of significant digits to facilitate the presentation. 



III. LOW ENERGY SPECTRA AND 
CORRELATION FUNCTIONS 

The ground state energies per spin along with their 
degeneracies are listed in Table HI They are plotted as 
a function of n in Fig. [21 The bigger clusters achieve 
the lowest energy, while C26 and C28 the highest. Prox- 
imity of polygons of the same kind is crucial for energy 
minimization [?|. C26 and C28 do not have hexagons ad- 
jacent to each other (Fig. [TJ, and their energy per spin 
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is even higher than the corresponding, energy for the do- 
decahedron, which has no hexagons [8(. They show that 
increase of the number of hexagons does not necessar- 
ily lead to lower energy per spin. For the bigger clus- 
ters the hexagons approach each other more and more 
and lower the energy, with 632,77 having the lowest. Its 
structure is such that there are two groups of three adja- 
cent hexagons each (Fig. [TJf ) ) , and all hexagon- hexagon 
bond correlations are strong (Table [VIII|) . In addition, 
there is a set of three pentagon-pentagon bonds in the 
middle of the cluster that have a very strong singlet char- 
acter with a correlation value equal to -0.60548. Such a 
strong correlation is not achieved even within hexagons 
in any of the clusters, and points to the importance of 
the specific structure geometry rather than the symme- 
try for the minimization of the energy, at least for the 
fullerene clusters with small n considered in this paper. 
The two hexagon groups and the strong singlet-like pen- 
tagon bonds are correlated weakly with the rest of the 
spins. 

The low energy spectra are shown in Tables Hill VIII and 
Fig. [3] Except from 62s, the ground and first excited 
states are singlets, with the ground states non-degenerate 
and the first excited states doubly degenerate except from 
C24 where it is non-degenerate. For 630 and 632,7 the 
second excited state is also a singly degenerate singlet. 
As seen in Figs. |3][d) and[3£e) the low energy spectra of 
C30 and 632,7 are similar. There is a non-degenerate sin- 
glet followed by two closely spaced singlets with the same 
degeneracy. The two clusters belong to different symme- 
try groups, but they both have their hexagons forming a 
band in the middle. Spin inversion symmetry is however 
opposite for the six lowest states in the spectra (Tables 
IVl and rVT|) . C26 and 6*32,77 on the other hand have the 
same spatial symmetry, D^h . In their low energy spectra 
(Tables QTI] and IVII[) the ground states belong to differ- 
ent irreducible representations, and their spin inversion 
symmetry is also different. The first excited states be- 
long to the same irreducible representation, however the 
spin inversion symmetry is still different. This is in con- 
trast to the smallest dodecahedral Ih symmetry clusters, 
the icosahedron (not of the fullerene type but made only 
of triangles) and the dodecahedron, which have similar 
structure and relative spacing of the levels in their low en- 
ergy spectra even though they comprise of different poly- 
gons, therefore symmetry is a strong determining factor 
for their properties. This result points to the conclusion 
that spatial symmetry is not the only factor determining 
the low energy properties in general for the fullerenes. 

C28 differs from the rest of the clusters in that its 
ground state is a doubly degenerate singlet @ , while the 
first excited state is a closely spaced triply degenerate 
triplet (Table HV|) . Then triply and doubly degenerate 
singlets follow, and the first non-degenerate state which 
has S = 2. In no other cluster an S = 2 state lies so 
low in the excitation spectrum. The ground state dou- 
blet belongs to the E s representation, which transforms 
as the pair (x 2 — y 2 , 2z 2 — x 2 — y 2 ) of Cartesian tensors 



. In contrast to [9( , here we find the correlations to be 
the same for both degenerate ground states. The ground 
state nearest neighbor correlation functions are shown in 
Table IVIIII for all the clusters. Except from the intra- 
pentagon correlation of 6*32,77 same hexagon correlations 
are the strongest. 

The nearest-neighbor correlation functions for the low- 
est singlet excitations are listed in table IIXI The two 
smallest clusters and 6*32,77 lower the energy of the 
hexagon-pentagon bonds while increasing the energy of 
the rest (only one of the same-hexagon bonds of 626 
lowers very weakly). On the contrary, C30 behaves the 
other way, except from the intra-hexagon correlation in 
the middle of the cluster that refers to a common side 
of two hexagons, which weakens. Its two singlet excited 
states are very close in energy and in the behavior of the 
correlation functions, even though they are of different 
multiplicity, which is also true for 632,7. The latter al- 
ters its same-hexagon bonds to generate the two excited 
singlets, even though the strongest one that refers to a 
common side of two hexagons does not change signifi- 
cantly. 

The nearest-neighbor correlation functions for the first 
triplet excitation are shown in Table [X] The values for 
632,77 change relatively little compared to the ground 
state, and the pentagon-pentagon correlation is getting 
even stronger. Similarly, 630 shows little change except 
from the hexagon-pentagon correlation that gets weaker. 
For 624 the triplet is mostly due to the decrease of the 
intra-hexagon correlation function, while in 626 the inter- 
hexagon correlation is getting stronger. In 62s only the 
inter-hexagon correlation < Si ■ S5 > changes relatively 
weakly compared to the other correlations, and correla- 
tion < Sio ■ Sii > is particularly weak. Finally, for 632,7 
there are strong changes for all correlations except from 
< 1S5 • <Si6 >■ 



IV. GROUND STATE MAGNETIZATION 

The ground state magnetization as a function of an 
external field has typically a step-like structure, with a 
AS = I discontinuity at fields where the ground state 
switches between adjacent S sectors. Frustration though 
can lead to magnetization discontinuities with AS > 1, 
where a particular S sector never necomes the ground 
state in a field. Such is the case for the icosahedral sym- 
metry J/j clusters, where the sector S — § — 5 with five 
flipped spins from saturation never includes the ground 
state [7] . The number of discontinuities is more than one 
at the classical level S — > 00, and also for the dodecahe- 
dron (n = 20) and Si — 1 where the calculation of the 
lowest energy state is computationally feasible for all S 
sectors (for Ih clusters with n > 20 the lowest energy 
state calculation is only possible for very high S even for 
Si = |). 

The lowest states for all the S sectors along with their 
degeneracies and the irreducible representation to which 
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they belong are listed in Table IXII (the saturation fields 
are listed in Table [I]). The corresponding reduced magne- 
tization M = — curves are shown in Fig. HI Unlike the 

nsi ° 1 — 1 

icosahedral symmetry case, no discontinuities are found. 

For some values of M there are plateaux, where a par- 
ticular 5 sector contains the ground state for a wider 
range of fields than the neighboring sectors. The most 
pronounced appears for C29, and 5 = 2, where M — 
0.14286 (Fig. He))). 

It is again hard to draw correlations between symme- 
try and the response in a magnetic field. In Figs. BJa) 
and [5ja) there is a correlation of the plateau-like fea- 
tures of the magnetization curve of C24 with stronger 
values of the intra-hexagon bonds. For C28, where the 
singlet-triplet gap is very small, there are stronger intra- 
hexagon correlation functions for the few low-lying 5 > 
sectors compared to the singlet case (Fig. EJc)). There 
are 5 sectors that contain the ground state for a very 
narrow range of the field (Fig. HJc)), and looking at Ta- 
ble IXII their lowest state belongs to three-dimensional 
irreducible representations. For C 30 , sectors 5 = 3 and 

4 have very strong same-hexagon correlations, stronger 
than the ones in the ground state (Fig. Eld)). For 6*32,7 
there are strong same-hexagon correlations for the low-5 
sectors, but for higher 5 intra-pentagon correlations are 
the strongest (Fig. Etc)). In both cases, the strength of 
these correlations does not change significantly with the 

5 value. 5 = 12 is the first sector that restores the same- 
hexagon correlations as the strongest, and it is the ground 
state for a narrow field window (Fig. [He)). In the case of 
C32.11 there are low-5 sectors where the intra-pentagon 
correlation is very strong (Fig. EJf)). For 5 — 1 and 3 it 
is even stronger than the 5 = value. For C26 the 5 = 4 
and 6 sectors are ground states for a narrow range of the 
field (Fig. 0Jb)) and same-hexagon correlations are weak 
(Fig. E^b)). Finally, for C24 the single spin-flip subspace 
has the ground state for a very narrow window of the 
field (Fig. (Ha)), with the two spin-flip subspace having 
a plateau and the strongest intra-hexagon correlations 
relative to its neighboring 5 sectors (Fig. Et a ))- 

For C26 and C28 the sector with a single spin flip from 
saturation is degenerate (Table IXip . However, the spin 
flips are not confined on the hexagons except from the 
singly degenerate state of C26, therefore there is no ana- 
log with the high magnetization localized magnon states 



discussed in Ref. [l2j . 



V. CONCLUSIONS 

The low energy spectrum and the magnetic response 
of the AHM have been calculated on a series of clusters 
with the connectivity of the fullerenes and a number of 
sites ranging from 24 to 32. Frustration and connectivity 
have a signature on the low energy spectrum with singlet 
ground and low energy excited states, the only excep- 
tion being the 28-site cluster where the ground state is 
a doubly degenerate singlet and the first excited state 
a triplet. Frustration is minimal when pentagons and 
hexagons minimize their interference in the clusters by 
being placed adjacent to polygons of the same kind. The 
magnetization as a function of an external field exhibits 
plateaux features, the most pronounced for n = 28 and 
5 = 2. Unlike the icosahedral Ih symmetry clusters [7|,|8j], 
spatial symmetry is not the sole determining factor of the 
magnetic properties of the clusters and their connectivity 
appears to be important as well. It is desirable to investi- 
gate clusters of higher n to gain more insight on the cor- 
relation between symmetry, connectivity and magnetic 
properties, however this is very challenging with present 
day computational means, at least for the low energy 
properties. For low n the competition between unfrus- 
trated hexagons and frustrated pentagons is strong, but 
even for high n the pentagon influence can be important, 
as in the case of Ih symmetry [7J • 
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TABLE I: Symmetry and ground state properties for the five 
clusters. — 1 is the ground state energy per spin and mult, is 
the state's multiplicity. h sat is the saturation field. 



sites 
n 


symmetry 
group 


number of 
symmetry 
operations 


Bo 

n 


mult. 


hsat 


24 


D 6d 


24 


-0.48831 


1 


4 + 72 


26 


D 3h 


12 


-0.48496 


1 


4 + ^/2 


28 


T d 


24 


-0.48482 


2 


4 + V2 


30 


D 5h 


20 


-0.49625 


1 


3 + s/7 


32 (I) 


D 3d 


12 


-0.49597 


1 


9+V5 


32 (II) 


D zh 


12 


-0.49804 


1 


5.61050 



TABLE IV: Low energy spectrum for C 2 g. Notation as in 
table H] 



E 


mult. 


irrep. 


s 


E 


mult. 


irrep. 


S 


-13.57486 


2 


E s 


o 


-13.22866 


3 


T 2 , a 


1 


-13.55978 


3 


T 2 ,a 


1 


-13.21762 


1 


A\ a 





-13.50468 


3 


T 2 , s 





-13.19914 


3 


T 2 , s 





-13.49099 


2 


E s 





-13.19486 


3 


Ti, 


1 


-13.42327 


1 


Ai, s 


2 


-13.18471 


2 


E 3 





-13.38652 


2 


E a 


1 


-13.17446 


3 


T 2>a 


1 


-13.36105 


3 


T 2 ,a 


1 


-13.14963 


1 


A 2 , s 





-13.29319 


3 


Ti,s 





-13.14309 


1 


Al,a 


1 


-13.27818 


3 


Tl, a 


1 


-13.11265 


3 




2 



TABLE II: Low energy spectrum for C24. E is the energy, 
mult, stands for the multiplicity of the state and irrep. for 
irreducible representation. S is the total spin, with each S 
state corresponding to 2S + 1 states with different projection 
of the total spin along the z axis S z . The spatial irreducible 
representation notation follows Ref. [nj. Indices s and a 
indicate the behavior under spin inversion, where s stands for 
symmetric and a for antisymmetric. A comma is introduced 
when necessary to avoid confusion between the notation for 
the spatial irreducible representation and the behavior under 
spin inversion. 



E 


mult. 


irrep. 


S 


E 


mult. 


irrep. 


S 


-11.71937 


1 


Bi >s 





-11.23652 


2 


E 2 , a 


1 


-11.70478 


1 


Ai tS 





-11.23043 


1 


A 1>s 





-11.46814 


2 


E3, a 


1 


-11.22852 


2 


El. a 


1 


-11.37244 


1 


A 2 , a 


1 


-11.17883 


2 


E'J, t a 


1 


-11.29779 


1 


B h s 





-11.16437 


2 


E 4 ,s 





-11.29326 


2 


E 2>s 





-11.16257 


2 


E\ >s 





-11.27888 


1 


B 2-a 


1 


-11.15826 


2 


Es,s 





-11.26341 


2 


Es^a 


1 


-11.13896 


2 


Ei.a 


1 


-11.24565 


2 


E5 t s 





-11.05126 


1 


A 1>s 






TABLE III: Low energy spectrum for C 2 a- Notation as in 
table HI 



E 


mult. 


irrep. 


S 


E 


mult. 


irrep. 


S 


-12.60898 


1 


A 2 . a 





-12.31502 


1 


A 2 , a 





-12.55739 


2 


E a 





-12.30577 


2 


E s 


1 


-12.49297 


2 


E s 


1 


-12.28884 


2 


E s 


1 


-12.46174 


1 


A2, a 





-12.26375 


2 


E a 





-12.44862 


1 


a'L 





-12.21057 


1 


A\. s 


1 


-12.42682 


2 


E a 





-12.16192 


2 


E s 


1 


-12.42129 


2 


E s 


1 


-12.15183 


2 


A\. a 





-12.38217 


1 


A±.s 


1 


-12.14598 


2 


Es 


1 


-12.37889 


1 


A\. s 


1 


-12.13819 


1 


A 2 , a 


2 



TABLE V: Low energy spectrum for C30. Notation as in table 

im 



E 


mult. 


irrep. 


s 


E 


mult. 


irrep. 


S 


-14.88742 


1 


A 2 . a 





-14.47476 


1 


A 2 , a 





-14.83815 


2 


n 

E 2 , a 





-14.46287 


2 


E 2 ,s 


1 


-14.82517 


1 


A 2 , a 





-14.45130 


2 


Ei,s 


1 


-14.62495 


1 


A ha 


1 


-14.33694 


2 


1 

El, a 





-14.60458 


2 


E 2 , 3 


1 


-14.27711 


2 


E\, a 





-14.59928 


1 


tt 

A hs 


1 


-14.26727 


1 


A 2 , a 





-14.51907 


2 


Ei, s 


1 


-14.23932 


1 




1 


-14.50190 


2 


El, a 





-14.17268 


2 


E 2 , a 


1 


-14.48316 


2 


E 2>a 


1 


-14.16221 


1 


A 2>a 


2 



TABLE VI: Low energy spectrum for C32,/. Notation as in 
table HI 



E 


mult. 


irrep . 


S 


E 


mult . 


irrep . 


s 


-15.87092 


1 







-15.49168 


1 




1 


-15.81199 


2 


p 





-15.47080 


1 







-15.80648 


1 







-15.46339 


1 


Al a 


1 


-15.67299 


2 


-Eu,a 


1 


-15.45531 


2 




1 


-15.59634 


2 


Eg.a 


1 


-15.44665 


2 







-15.57987 


1 


A2u,a 


1 


-15.44513 


2 




1 


-15.52875 


2 


Eg,a 





-15.43356 


1 


A2g,a 


1 


-15.51890 


1 


Alg, S 





-15.38287 


1 


Alg,B 


2 


-15.49978 


2 


E u ^ a 


1 


-15.38034 


1 


A2u,a 


1 



TABLE VII: Low energy spectrum for Ci2,n- Notation as in 
table El 



E 


mult. 


irrep. 


S 


E 


mult. 


irrep. 


S 


-15.93723 


1 


A\ tS 





-15.50045 


1 


A' 2 \ a 


1 


-15.81192 


2 


E s 





-15.49288 


1 


A\ t s 





-15.77366 


1 


A' 2 ,a 


1 


-15.46219 


1 


A\ :S 





-15.73368 


1 


A\^ s 





-15.45437 


1 


A' 2 ,a 


1 


-15.63730 


2 


E a 


1 


-15.45317 


1 


^4-1, a 


1 


-15.60167 


2 


E a 


1 


-15.42185 


2 


E s 





-15.57485 


2 


E s 





-15.39363 


1 


Ala 


1 


-15.56589 


2 


E a 


1 


-15.38231 


1 


A'la 


1 


-15.54070 


2 


E a 


1 


-15.38020 


2 


E a 


1 



TABLE VIII: Distinct nearest-neighbor correlation functions 
for the ground states. 



sites n 


intra-hcxagon 


inter-hexagon 


hexagon-pentagon 


intra-pentagon 


24 


< Si- S 2 >= -0.40325 




<Si- S 7 >= -0.20285 


< S 7 ■ S s >= -0.37051 


26 


< S 5 • S 7 >= -0.33858 

< 5 5 • Su >= -0.42436 


< 5n • 5i2 >= -0.10317 


< S 2 ■ 5 5 >= -0.26523 


< Si • 5 2 >= -0.33215 


28 


< Si- S 2 > - -0.35418 


< Si • 5 5 >= -0.23883 


< S 2 • 5*3 >= -0.30347 




30 


< Se ■ S 7 >= -0.36221 

< S 7 ■ Sis >= -0.35466 




<Si- S 6 >= -0.24081 


< Si • S 2 >= -0.34618 


32 (D 3d ) 


< S$- S 6 >= -0.37205 

< S 5 ■ Sie >= -0.32629 

< S u ■ S12 >= -0.24979 

< Su ■ 5i9 >= -0.45409 




< 5 2 • S 6 >= -0.30381 


< Si • S 2 >= -0.30903 


32 (D 3h ) 


< Si- S 2 >= -0.36273 

< S 2 - S 5 >= -0.35602 

< 5 5 ■ 5n >= -0.41480 


< 5n • S22 >= -0.24352 


< 5 5 • 5 M >= -0.16367 


< S14 ■ S15 >= -0.60548 



TABLE IX: Distinct nearest-neighbor correlation functions 
for the singlet excited states within the singlet-triplet gap. 



sites n 


intra-hcxagon 


inter-hexagon 


hexagon- p ent agon 


intra-pentagon 


24 


< Si ■ 5 2 >= -0.37785 




< 5i ■ 5 7 >= -0.25586 


< Sr ■ 5 8 >= -0.34169 


26 


< S 5 ■ SV >= -0.28615 

< 5 5 • Su >= -0.42640 


< 5n • 5i2 >= -0.084111 


< 5 2 • 5 5 >= -0.30251 


< 5i • 5 2 >= -0.30687 


30 


< 5 6 ■ S 7 >= -0.37609 

< 5*7 • 5i 8 >= -0.33026 




< 5i • 5 6 >= -0.21102 


< 5i • 5 2 >= -0.35548 


30 


< 5 6 • 5 7 >= -0.37784 

< 5*7 • 5i 8 >= -0.32507 




< 5i • 5 6 >= -0.20832 


< 5i • 5*2 >= -0.35598 


32 (D 3d ) 


< 5 5 • 5 6 >= -0.39744 

< 5*5 ■ 5i6 >= -0.29843 

< 5n • 5i2 >= -0.27220 

< 5n • 5i9 >= -0.45947 




< 5 2 • 5 6 >= -0.30566 


< Si ■ 5 2 >= -0.29805 


32 (D 3d ) 


< 5 5 • 5 6 >= -0.40407 

< 5*5 • 5i6 >= -0.30669 

< 5n ■ 5i2 >= -0.27231 

< 5n • 5i9 >= -0.45357 




< 5*2 • 5*6 >= -0.29314 


< 5i • 5*2 >= -0.30482 


32 (D 3h ) 


< 5i • 5* 2 >= -0.36028 

< 5* 2 • 5*5 >= -0.34068 

< 5*5 • 5n >= -0.41156 


< 5n • 5*22 >= -0.16865 


< 5 5 • 5i4 >= -0.24184 


< 5i4 • 5i5 >= -0.40508 



TABLE X: Distinct nearest-neighbor correlation functions for 
the first triplet excited states. 



sites n 


intra-hcxagon 


intcr-hexagon 


hexagon- p ent agon 


intra-pentagon 


24 


< 5i • 5 2 >= -0.38523 




< 5i ■ 5 7 >= -0.20358 


< 5 7 ■ 5g >= -0.36686 


26 


< 5 5 • 5 7 >= -0.34294 

< 5 5 • 5ii >= -0.40945 


< 5n • 5i2 >= -0.13733 


< 5 2 • 5 5 >= -0.25900 


< 5i • 5 2 >= -0.33365 


28 


< 5i • 5 2 >= -0.32073 

< 5*4 • 5*5 >= -0.42032 

< 5*2 • 5io >= -0.42541 


< 5i • 5 5 >= -0.22540 
< 5io ■ 5ii >= -0.010665 


< 5*2 • 5*3 >= -0.33309 

< 5*3 ■ 5*4 >= -0.16012 




30 


< Sa ■ S 7 >= -0.35985 

< 5*7 • Sis >= -0.34801 




< 5i • 5 6 >= -0.22022 


< Si ■ 5 2 >= -0.34857 


32 (D 3d ) 


< 5 5 • 5 6 >= -0.40953 

< 5*5 • 5i6 >= -0.33089 

< 5n • 5i2 >= -0.28203 

< 5n • 5i9 >= -0.42064 




< 5 2 • 5 6 >= -0.25415 


< 5i • 5 2 >= -0.32989 


32 (D 3h ) 


< 5i • 5 2 >= -0.35863 

< 5 2 • 5* 5 >= -0.35762 

< 5*5 ■ 5ii >= -0.40890 


< 5n • 5 22 >= -0.24266 


< 5 5 • 5 W >= -0.15494 


< 5 W • 5is >= -0.61209 



TABLE XI: Lowest energies (-Eo), multiplicities (mult.) and 
corresponding irreducible representations (irrep.) in each to- 
tal spin S sector for the clusters. The spatial irreducible 
representation notation follows Ref. [IlJ. Indices s and a 
indicate the behavior under spin inversion, where s stands 
for symmetric and a for antisymmetric. It is only possible 
to calculate them for states lying relatively low or high in 
the energy spectrum. A comma is introduced when necessary 
to avoid confusion between the notation for the spatial irre- 
ducible representation and the behavior under spin inversion. 







C24 


C24 


C26 


C26 


C26 


C28 


C28 


C*28 


s 


E 


mult. 


irrep. 


E 


mult. 


irrep. 


Eo 


mult. 


irrep. 





-11.71937 


1 




-12.60898 


1 


A 2 ,a 


-13.57486 


2 


E s 


1 


-11.46814 


2 




-12.49297 


2 


E s 


-13.55978 


3 


T2,a 


2 


-10.92259 


1 




-12.13819 


1 


A 2 ,a 


-13.42327 


1 


Ai tS 


3 


-10.13390 


1 


B2,a 


-11.39120 


1 


Als 


-12.45085 


3 


T2,a 


4 


-8.88178 


1 


Ai j3 


-10.16420 


1 


A 2 ,a 


-11.45102 


1 


Ai tS 


5 


-7.28948 


2 


E 5 


-8.83042 


1 


A'ls 


-10.18889 


1 


Al, a 


6 


-5.49575 


1 


A 1 


-7.01596 


1 


a' 2 


-8.51700 


1 


Ai 


7 


-3.46016 


1 


B 2 


-5.07125 


1 


a\ 


-6.56848 


3 


T 2 


8 


-1.24214 


1 


Si 


-2.91337 


2 


E 


-4.52626 


2 


E 


9 


1.10955 


1 


A 2 


-0.62134 


2 


E 


-2.32108 


1 


A 1 


10 


3.59067 


1 




1.75066 


1 


a' 2 


-0.049533 


1 


Ai 


11 


6.29289 


2 


Ez,a 


4.36198 


1 


A! 


2.48354 


3 


T 2 


12 


9 


1 


Ai tS 


7.04289 


1,2 


A 2 . a ,E 


5.11292 


2 


E 


13 








9.75 


1 


1 

A 1>s 


7.79289 


1,3 


Al, a ,T 2 


14 














10.5 


1 


Ai, s 






C30 


C30 


Cz2,I 


C32,/ 


C32.1 


C32,II 


C32, ji 


Cs2,II 


S 


Eo 


mult. 


irrep. 


Eo 


mult. 


irrep. 


Eo 


mult. 


irrep. 





-14.88742 


1 


n 

A 2 .a 


-15.87092 


1 




-15.93723 


1 


n 

Ai, a 


1 


-14.62495 


1 


A 1<3 


-15.67299 


2 


E ua 


-15.77366 


1 


A 2 .a 


2 


-14.16221 


1 


1 

A 2 .a 


-15.38287 


1 


Algs 


-15.35876 


1 


t 

Ai 


3 


-13.47738 


2 


E2,s 


-14.72334 


1 


A 2g 


-14.68665 


1 


A 2 


4 


-12.52542 


1 


A 


-13.83231 


1 


Alg 


-13.80765 


1 


A 


5 


-11.26275 


1 


4 


-12.66407 


1 


A 2g 


-12.59064 


1 


A'i 


6 


-9.71967 


2 


E 2 


-11.25111 


1 


Alg 


-11.12132 


1 


11 

Ai 


7 


-7.94218 


2 


E 2 


-9.56173 


1 


A 2g 


-9.42233 


1 


A'i 


8 


-6.07738 


1 


II 

A 2 


-7.65604 


1 


Ai g 


-7.61188 


1 


1 

Ai 


9 


-3.98134 


1 


A 1 


-5.59828 


1 


A 2g 


-5.58108 


1 


A 2 


10 


-1.69706 


1 


A 2 


-3.40688 


1 


Alg 


-3.43626 


1 




11 


0.66164 


1 


a\ 


-1.06488 


1 


A 2g 


-1.14953 


1 


A 2 


12 


3.08988 


1 


A 2 


1.37454 


1 


A 2g 


1.25805 


1 


A 


13 


5.68792 


1 


a\ 


3.84610 


1 


A 2g 


3.79378 


1 


A 2 


14 


8.42712 


1 


A'l 


6.44732 


1 


Ai g 


6.42049 


1 


A'i 


15 


11.25 


1 




9.19098 


1 


A 2g 


9.19475 


1 


A'2 


16 








12 


1 


Alg, a 


12 


1 


A\, 3 



9 




(c) (d) 




(e) (f) 



FIG. 1: Projection of the clusters on a plane: (a) C24, (b) 
C 2 6, (c) C 2 8, (d) C30 (e) C32,/ and (f) 632,11- For C30 there 
is also a top (bottom with the dashes) view. The black circles 
are spins s». The solid lines are antiferromagnetic bonds J. 
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FIG. 2: Ground state energy per spin — a as a function of the 
number of spins n. The value for n — 20 is taken from 
The lowest energy for n = 32 is for cluster C32JI, and the 
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FIG. 3: Low energy E spectrum of total spin S states and its 
multiplicity: (a) C24, (b) C 2 6, (c) C28, (d) C30 (e) C32,/ and 
(f) C$2,11- Solid lines: S = 0, dashed lines: S = 1, dotted 
lines: 5 = 2. 
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FIG. 4: Reduced ground state magnetization M = as 
a function of magnetic field h: (a) C24, (b) C26, (c) C28, (d) 
C30, (e) C32,/ and (f) C32,//. M is the total spin 5 normalized 
to the number of sites n and the magnitude of spin s;. M has 
no units and h is in units of energy. 
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FIG. 5: Distinct correlation functions for the lowest energy 
state in each total spin S sector: (a) C24, (b) C26, (c) C28, 
(d) C30, (e) C32,/ and (f) C32,//. The reduced magnetization 
M = is the total spin S normalized to the number of 

sites n and the magnitude of spin Sj. < Si ■ Sj > is in units 
of energy and M has no units. < Si ■ Sj >: A, <, >, V: 
intra-hexagon, o, o: inter-hexagon, +, □: hexagon-pentagon, 
x: intra-pentagon. 



